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Abstract. Massive and intermediate mass stars play a crucial role in astrophysics. Indeed, massive stars are 
the main producers of heavy elements, explode in supernovae at the end of their short lifetimes, and may be 
the progenitors of gamma ray bursts. Intermediate mass stars, although not destined to explode in supemovae, 
display similar phenomena, are much more numerous, and have some of the richest pulsation spectra. A key to 
understanding these stars is understanding the effects of rapid rotation on their structure and evolution. These ef¬ 
fects include centrifugal deformation and gravity darkening which can be observed immediately, and long terms 
effects such as rotational mixing due to shear turbulence, which prolong stellar lifetime, modify chemical yields, 
and impact the stellar remnant at the end of their lifetime. In order to understand these effects, a number of 
models have been and are being developed over the past few years. These models lead to increasingly sophisti¬ 
cated predictions which need to be tested through observations. A particularly promising source of constraints is 
seismic observations as these may potentially lead to detailed information on their internal structure. However, 
before extracting such information, a number of theoretical and observational hurdles need to be overcome, not 
least of which is mode identification. The present proceedings describe recent progress in modelling these stars 
and show how an improved understanding of their pulsations, namely frequency patterns, mode visibilities, line 
profile variations, and mode excitation, may help with deciphering seismic observations. 


1 Introduction 

Massive stars play an important role in astrophysics. In¬ 
deed, they are the main producers of heavy elements, dom¬ 
inate the chemical evolution of galaxies, and end their lives 
in supernovae. Although intermediate mass stars do not ex¬ 
plode in supernovae, they present many of the same phe¬ 
nomena as massive stars, and furthermore can have very 
rich pulsation spectra. Therefore, understanding the struc¬ 
ture and evolution of these stars, in particular their chem¬ 
ical evolution and internal transport processes, is crucial. 
Unlike their less massive counterparts, these stars hardly 
have any surface convection (apart from possibly a thin 
near-surface convective layer). Consequently, they are not 
spun down and tend to be rapid rotators {e.g. ||T|). Hence, 
in order to understand these these stars, one needs to un¬ 
derstand the impact of rapid rotation on stellar structure, 
chemical evolution, and pulsations. 

2 Impact of rotation on stellar structure 
and evolution 

Rapid stellar rotation has both immediate effects on stellar 
structure as well as long range effects on stellar evolution. 

2.1 Structural changes 

One of the first and most obvious effects of rapid rotation 
is stellar deformation. Indeed, the centrifugal force causes 
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the equator to bulge, as confirmed by recent interferomet¬ 
ric observations (e.g. fJl, |0, ||4|). This effect is important 
because it means these stars can no longer be modelled us¬ 
ing a ID approach but require a 2D approach. 

A second effect, which also has immediate consequences, 
is gravity darkening. Indeed, in rapidly rotating stars, the 
equator is less luminous than the poles, due to a smaller 
vertical temperature gradient. This effect has also been ob¬ 
served through interferometry (e.g. 0, 0, Q) and is im¬ 
portant because it causes the position of these stars in an 
HR diagram to depend on their inclination (e.g. 0, ii). 

Various approaches have been used to model this ef¬ 
fect, the first being von Zeipel’s law IfTOl . This approach 
corresponds to the simple power law Teff oc where 
P - 0.25. The value p - 0.8 has also been proposed for 
stars with a convective envelope IfTTI . More recently, Es¬ 
pinosa Lara and Rieutord came up with a new law which 
is based on two assumptions: first, that the luminous flux is 
parallel to the effective gravity, and second that the diver¬ 
gence of the luminous flux is zero (which corresponds to 
assuming no production of energy in the envelope). This 
leads to an analytical solution which does not depend on 
any free parameters, and which compares very favourably 
with full 2D numerical simulations from the ESTER code 

Ha. 


2.2 Baroclinic effects 

As opposed to their non-rotating counterparts, rapidly ro¬ 
tating stars have a baroclinic structure. In other words, lines 
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of constant pressure, temperature, and density no longer 
coincide. This produces differential rotation which in turn 
leads to meridional circulation due to viscosity US). Cur¬ 
rently, only the ESTER code is able to self-consistently 
calculate such flows in realistic 2D stellar models (M, 

my 

As a result of various instabilities, baroclinic flows will 
then cause turbulence and enhance transport processes, which 
modify the stellar lifetime and chemical yields {e.g. my 
Predictions from ID stellar evolution codes including a 
shellular rotation profile based on El are in better agree¬ 
ment with various observations, such as the chemical en¬ 
richment of OBA stars, the number ratio of red to blue 
supergiants etc. (see El and references therein). Recent 
observations of N enrichment in O and B stars ( lfT9l . ij^Ol l 
have, however, shown discrepancies with theoretical ex¬ 
pectations. Indeed, there is a lack of correlation between 
N enrichment and equatorial velocity, whereas rotational 
mixing is needed in theoretical models to produce N en¬ 
richment. It is not clear what could cause this, although 
various solutions have been proposed, such as transport of 
chemical species by waves CTI . 

2.3 Impact on convection zones 

Rapid rotation also affects convection zones. Eor instance, 
1221 produced a rapidly rotating model with a convectively 
unstable equatorial belt. A year later, ||23l, using a ID for¬ 
malism, found that rotation favours convection in stellar 
envelopes, especially at the equator. Nonetheless, it remains 
an open question as to what prescription for convection 
should be used in a 2D rapidly rotating model, and is cur¬ 
rently the main reason why the ESTER code is unable to 
model low mass stars. However, ongoing work by Il24l 
may provide an answer to this question. Indeed, they are 
developing a new convection code based on unstructured 
meshes, capable of handling centrifugal deformation. Ini¬ 
tial comparisons with benchmarks from the ASH code show 
promising results. 

2.4 Summary 

In summary, rotation causes many new phenomena which 
affect stellar structure, transport processes, mixing, and evo¬ 
lution. Although much progress has been made in our the¬ 
oretical modelling of these processes, there remains large 
uncertainties due to their complexity and to discrepancies 
with current observations. As a result, it is necessary to 
constrain these processes through further observations, and 
asteroseismology is curi'ently the best way of doing this as 
it is the only way we have to probe the internal structure of 
stars. However, before being able to carry out asteroseis¬ 
mology, one needs to be able to model the effects of rapid 
rotation on stellar pulsations, which, as will be described 
in the next section, is by no means trivial. 


3 Impact of rotation on stellar pulsations 

In the non-rotating case and in the absence of phenomena 
which break spherical symmetry, such as magnetic fields, 
stellar pulsation modes with the same radial order, n, and 


harmonic degree, and 1, but differing azimuthal orders, m, 
have the same frequencies and are thus degenerate. How¬ 
ever, it has been known for many years that rotation lifts 
this degeneracy ||25|, in much the same way that the Zee- 
man effect splits absorption lines. The first order effect of 
rotation can be split into two terms. The first is simply an 
advection term which leads to a sort of Doppler shift of 
the frequencies on non-axisymmetric modes. The second 
comes from the effects of the Coriolis force and is char¬ 
acterised by the Ledoux constant. Inversions of solar or 
stellar rotation profiles rely on this first order approxima¬ 
tion (e.g. m, El, El). At more rapid rotation rates, 
one needs to include higher order effects of rotation (e.g. 
EH, ESI), ultimately apply a 2D approach (e.g. ED, 
ED). In what follows, we will describe the effects of rapid 
rotation on stellar pulsations, as based on the latest 2D nu¬ 
merical simulations. 


3.1 Gravito-inertial modes 

At the lower end of the frequency spectrum in non-rotating 
stars, there are gravity (or g-) modes for which the restor¬ 
ing force is buoyancy. When the star is rotating, buoy¬ 
ancy combines with the Coriolis force, thereby leading to 
gravito-inertial modes. The type of stars which display gravito- 
inertial modes are y Dor stars, SPBs, and Be stars. If the 
mode owes its existence to the Coriolis force, then it is 
known as an inertial mode. An extensive literature on iner¬ 
tial modes, and singular gravito-inertial modes exist (e.g. 

ED, ED, ED, ED, ED), but in what follows, we will 

focus on modes which become classical g-modes in the 
non-rotating limit. 

According to TassouTs asymptotic formula ED, grav¬ 
ity modes in non-rotating stars are evenly spaced in period, 
for a given 1 value. When the star is rotating, the period 
spacing also depends on m and on the ratio V - ^ which 
characterises the effects of the Coriolis force on the pul¬ 
sation mode, where Q is the rotation rate and u) the pul¬ 
sation frequency. Such a relation was first established us¬ 
ing the traditional approximation (EH, iHOl l and has been 
subsequently confirmed through full 2D calculations, al¬ 
though the centrifugal deformation slightly affects some of 
the modes BTIl . 

In terms of mode geometry, it is extremely important to 
distinguish between sub-inertial modes, for which w < 2Q, 
and super-inertial modes, where w > 2Q. Indeed, in the 
sub-inertial regime, modes are confined to an equatorial 
region thanks to critical surfaces, as described in Il42l and 
ED- In the super-inertial regime, gravito-inertial modes 
tend to behave like their non-rotating counterparts, except 
for one notable exception - rosette modes. Such modes 
where discovered by BTI using a 2D numerical approach, 
and have subsequently been studied in detail by Il43l . Il44l . 
ED, ED using a variety of techniques. 

3.2 Acoustic modes 

At the upper end of the frequency spectrum, there are acous¬ 
tic (or p-) modes for which the restoring force is pressure. 
Given that these modes tend to be located in the outer 
portions of the star, they are most affected by centrifugal 
deformation, which is greatest in such regions. One way 
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of characterising the effects of the centrifugal deforma¬ 
tion is by taking the ratio of the characteristic length scale 
of the deformation, which is roughly proportional to Q^, 
to the mode’s wavelength, which is proportional to w 
Hence, high frequency modes are more affected than low 
frequency modes, due to their smaller wavelength, d Scuti 
and p Cephei stars oscillate in the acoustic domain and typ¬ 
ically tend to be rapid rotators. 

At rapid rotation rates, the acoustic modes subdivide 
into different classes of modes, each with their own charac¬ 
teristic geometry and independent frequency organisation. 
This was first shown in poly tropic models by If47l . Il48l 
using 2D numerical simulations and ray dynamics and was 
then extended to more realistic models using 2D numerical 
simulations by ESI. 

Among these different classes, island modes are the 
most important, since they are the most visible of the reg¬ 
ular (non-chaotic) modes. These modes are the rotating 
counterpart to modes with few nodes from one pole to 
another, i.e. with a small 1 - \m\ value. Their geometric 
structure closely follows a 2-period ray orbit and is char¬ 
acterised by new quantum numbers as illustrated in Fig. [1] 
Using these quantum numbers, it is possible to obtain a 
new asymptotic formula which describes the frequencies 
of these modes, the coefficients of which are related to 
travel time integrals (see PTTlI . Il^ , 1501 and ISTl '). Of par¬ 
ticular interest is the frequency spacing, Af,, between modes 
with consecutive h values. This spacing is, in fact, half the 
large separation, and roughly scales with the square-root 
of the stellar mean density, even a rapid rotation rates l52l . 



M = 25.0Mo r| = 0.6 a = 0.0 

cj = 307.7/xHz m=l 

Fig. 1. Island mode in a model based on the Self-Consistent Field 
method 1531. New quantum numbers are indicated hy the arrows. 


3.3 Mixed modes 

Mixed modes combine the characteristics of acoustic and 
gravity modes, and occur when there are avoided cross¬ 
ings between these modes. Such modes typically occur in 


evolved stars, but can also be found in unevolved rapidly 
rotating stars. Indeed, gravito-inertial modes are hardly af¬ 
fected by centrifugal deformation and hence tend to scale 

as where Rpoi is the polar radius. In contrast, 

acoustic modes and the characteristic spacing zlj? are roughly 
proportional to yfGMjV, where V is the volume of the 
star. As a result, the frequency domain of acoustic modes 
decreases compared with that of gravito-inertial modes, 
thereby causing the two to overlap. This favours avoided 
crossings and hence mixed modes, especially given that 
rotation allows coupling between modes with different har¬ 
monic degrees. 

In evolved stars, rotation may play an important role 
in mixed modes. Indeed, as shown by l54ll . rotation affects 
different members of a given rotation multiplet differently, 
even at relatively small rotation rates. This is because the 
avoided crossings which produce the mixed modes occur 
at slightly different values of the rotation rate for different 
members of the multiplet. As a result, the ratio of the p- and 
g-mode contributions will be different, thereby leading to 
different mode inertias, and the frequencies will no longer 
be evenly spaced. 


4 Asteroseismology 

Having described the effects of rapid rotation on stellar 
pulsations, we now turn our attention to the asteroseismic 
inferences which can be drawn from observed stellar pul¬ 
sations. In what follows, we will distinguish between aver¬ 
age or global asteroseismology, which focuses on the gen¬ 
eral characteristics of the pulsation spectrum, and detailed 
or “boutique” asteroseismology, which relies on identify¬ 
ing individual modes. 


4.1 Global asteroseismology 

As was done above, we first start with the lower end of 
the frequency spectrum. Recent observations have shown 
that the pulsation frequencies of gravito-inertial modes in a 
number of Be stars tend to clump together. As explained in 
various publications {e.g. El, El, El), the co-rotating 
frequencies of the pulsation modes are likely to be much 
smaller than the rotation rate. When viewed from an in¬ 
ertial frame, these frequencies are Doppler shifted by inQ 
as described above, thereby leading to separate clumps for 
each azimuthal order. Accordingly, the frequency differ¬ 
ences between these clumps provide an asteroseismic mea¬ 
sure of the rotation rate. However, a recent publication, 
ll58l . showed discrepancies between this seismic rotation 
rate and more classical measurements of Q based on spec¬ 
troscopy. 

At higher frequencies, recent studies of p-modes in rapidly 
rotating 6 Scuti stars have found recurrent frequency spac- 
ings using histograms of frequency differences or Fourier 
transforms of frequency spectra ( ll59l . IMI . IMl . ||62l, ll63l . 
lIMl l. Similar studies have also been carried out using theo¬ 
retical pulsation spectra based on realistic mode visibilities 
(El, El) have shown that in favourable cases, such 
recurrent spacings could correspond to the large frequency 
spacing or half its value {i.e. An), or to a multiple of the 
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rotation rate. Interpreting these spacings as the large sepa¬ 
ration, ll62ll then went on to constrain the mean density of 
a 6 Scuti star observed with CoRoT. 


4.2 Detailed asteroseismology 

In order to obtain tighter constraints on stellar properties, 
one needs to carry out detailed asteroseismology. However, 
up until now, it has proven to be very difficult to correctly 
match observed frequencies with theoretical modes, i.e. 
identify modes, in rapidly rotating stars due to the lack 
of simple frequency patterns, as is very well illustrated 
in Fig. 5 of ll67l . In order to overcome these difficulties, 
one can envisage various approaches. One can still try to 
look for the more complicated patterns through the use 
of Echelle diagrams ( ll62l . ll^ . Il68l l or by matching an 
asymptotic formula to the frequencies Il69l . The latter ap¬ 
proach will however be thrown off by the presence of chaotic 
modes which follow their own frequency distribution ll48l . 
Another possibility is to apply observational mode iden¬ 
tification techniques based on multi-colour photometry or 
spectroscopy. Such techniques have already been applied 
in slowly rotating stars, but they need to be adapted to more 
rapid rotation. 

In the photometric approach, one typically looks at the 
amplitude ratios or phase differences of a given pulsation 
mode in different photometric bands. The advantage of this 
approach is that the intrinsic mode amplitude factors out, 
thereby leaving information which is only based on mode 
geometry. As opposed to the non-rotating case, the ampli¬ 
tude ratios and phase differences depend both on the az¬ 
imuthal order and on the inclination {e.g. EQl, El]). How¬ 
ever, at a given inclination, similar amplitude ratios are ob¬ 
tained for modes with the same 1 and m but different radial 
orders Il72l . This is due to the similar surface geometric 
structure of these modes, as expected from asymptotic the¬ 
ory as]. Accordingly, by grouping together modes with 
similar amplitude ratios, one can find families of modes 
with similar quantum numbers, thereby constraining mode 
identification and characteristic spacings such as the large 
frequency separation or the rotation rate ll66] . 

The spectroscopic approach consists in obtaining a time 
series of high resolution spectra of a pulsating star to see 
how the absorption lines vary over time as a result of Doppler 
shifts from the oscillatory motions, and comparing this to 
theoretical predictions. Such an approach cannot currently 
be carried out from space and requires the use of ground 
based-telescopes, preferably in a network such as SONG. 
Previous and current spectroscopic observations of rapidly 
rotating pulsating stars include those described in Ea,E4), 
ES . Various studies have focused on obtaining theoretical 
predictions in a rapidly rotating context (EQ. EtI . ED) 
but more work is needed before these effects are incorpo¬ 
rated into mode identification tools such as FAMIAS 1791 . 

An important ingredient in predicting theoretical pho¬ 
tometric and spectroscopic signatures of pulsation modes 
is the variations of effective temperature. However, this 
quantity can only be calculated through non-adiabatic pul¬ 
sation calculations. Currently, the only studies that include 
these effects in rapidly rotating stars are those based on 
ifSOl , which use models based on a Chandrasekhar expan¬ 
sion and focus on low frequency modes, and ongoing work 
which involves using the TOP pulsation code ( 11311 . iH^l l 


with models from the ESTER code, an approach which is 
applicable to all modes. Figure |2] shows a 2D excitation 
map as well as some ID work integrals of stable and un¬ 
stable acoustic island-modes based on these latest calcu¬ 
lations. An added benefit of non-adiabatic calculations is 
knowing which modes are unstable, which could provide 
additional constraints on the mode identification. 




log(T) 


Fig. 2. (Upper panel) Two dimensional map of regions which 
excite (red/white) or damp (blue/black) a given island mode in an 
ESTER model. The radial coordinate in this plot is log T which 
stretches out the near surface regions where most of the excitation 
and damping takes places. (Lower panel) Work integrals for a 
multiplet of modes at £? = 0.4 As can be seen, some of the 
members of this multiplet are stabilised by rotation. 


5 Conclusion 

Rapid rotation plays a major role in the structure, evo¬ 
lution, and pulsations of massive and intermediate mass 
stars. Given that these stars are important for many do¬ 
mains of astrophysics, it is essential to understand the im¬ 
pact of rapid rotation. Recent theoretical developments are 
leading to more and more realistic models and a better un¬ 
derstanding of their pulsation modes. Nonetheless, many 

















The Space Photometry Revolution 


questions remain unanswered, especially in the light of 
new and spectacular observations. Progress has been made 
in interpreting the general characteristics of observed pul¬ 
sation spectra, but more work is needed in order to fully 
exploit the excellent data from current (MOST, CoRoT, 
Kepler, Brite) and future missions (TESS, PLATO). 
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